Displacement operator for quantum systems with position-dependent mass 



O 
(N 
» 

Q 

o 



Oh 
I 

^ . 



> 

00 

in 



R. N. Costa Filho/'B M. P. Almeida/ 'Q G. A. Farias/'l and J. S. Andrade Jr.^'l 

^ Departamento de Fisica, Universidade Federal do Ceard, 
Caixa Postal 6030, Campus do Pici, 60455-760 Fortaleza, Ceard, Brazil 

A translation operator is introduced to describe the quantum dynamics of a position-dependent 
mass particle in a null or constant potential. From this operator, we obtain a generalized form of 
the momentum operator as well as a unique commutation relation for x and p-,. Such a formalism 
naturally leads to a Schrodinger-like equation that is reminiscent of wave equations typically used to 
model electrons with position-dependent (effective) masses propagating through abrupt interfaces in 
semiconductor heterostructures. The distinctive features of our approach is demonstrated through 
analytical solutions calculated for particles under null and constant potentials like infinite wells in 
one and two dimensions and potential barriers. 
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The concept of non-commutative coordinates as a way 
of removing divergences in field theories through a mii- 
versal invariant length parameter was originally proposed 
by Heisenberg [jy]. This idea led to space-time quanti- 
zation 0, where a non-commutative space operator is 
defined to allow for the development of a theory invari- 
ant under Lorentz transformation, but not invariant un- 
der translations. Accordingly, the last property leads to 
continuous space-time coordinates. Since then, the need 
for a fundamental length scale became evident in differ- 
ent areas of physics, such as relativity string theory 
and quantum gravity Q. In the particular case of 
quantum mechanics, this minimum length scale yields a 
modification in the position momentum commutation re- 
lationship 

Previous studies have shown that the modification of 
canonical commutation relations or any modification in 
the underlying space typically results in a Schrodinger 
equation with a position-dependent mass [l^. This ap- 
proach has been rather effective in the description of elec- 
tronic properties of semiconductor ll| and quantum dots 
12| . Under this framework, mass is turned into an op- 
erator that does not commute with the momentum op- 
erator. This fact immediately raises the problem of the 
ordering of these operators in the kinetic energy operator 
[l3| . In this Communication, we introduce the concept 
of non-additive spatial displacement in the Hilbert space. 
This property not only changes the commutation relation 
for position and momentum, which leads to a modified 
uncertainty relation, but also reveals a Schrodinger-like 
differential equation that can be interpreted in terms of a 
particle with position-dependent mass, leading to a nat- 
ural derivation of the kinetic operator for such problems. 

Consider a well-localized state around x that can be 
changed to another well-localized state around x+a+jax 
with all the other physical properties unchanged, where 
the parameter 7 is the inverse of a characteristic length 
that determines the mixing between the displacement 
and the original position state. For 7 7^ 0, the displace- 
ment depends explicitly on the position of the system. 



while 7 = corresponds to a standard translation. This 
process can be mathematically expressed in terms of the 
operator Ty{a) as, 



T-y{a)\x) = \x + a + jax). 



(1) 



The composition of displacements through 7^ in terms of 
two successive infinitesimal translations results in 

T-,{dx')%{dx") = T-,{dx' + dx" -f jdx'dx"), (2) 

which clearly shows the non-additivity characteristic of 
the operator. It is also important to note that the inverse 
operator is given by. 



T-\dx)\x) 



dx 



1 + ^dx 



(3) 



Another relevant property is that 7^ becomes an identity 
operator when the infinitesimal translation goes to zero. 



lim Ty(dx) 

dx^O 



1. 



(4) 



At this point, we observe that the operator Tj{x) is the 
infinitesimal generator of the group represented by the 
so-called g-exponential function originally defined as , 



exp (a;) = [1 + (1 - (?)a 



(5) 



where expg(a) expg(6) — exp^[a + 6 -I- (1 — q)ab], and 
expo(a) = exp(a). The definition ([5]) represents a crucial 
ingredient in the mathematical formalism of the general- 
ized Tsallis thermostatistics flBI and its several applica- 
tions related with non-additive physical systems [16-18|. 
Associating Ty{dx) with the g-exponential, and expand- 
ing it to first order in dx leads to 



Tyidx) 



1 



ipjdx 



(6) 



where Pj is a generalized momentum operator and we are 
using the fact that momentum is a generator of transla- 
tion. Now, considering that 

xT-y{dx)\x) ^ (x + dx + jxdx)\x + dx + jxdx) , (7) 
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and 7''y{dx)x\x) = x\x -'r dx -\- 'yxdx), we obtain the com- 
mutation relation, 

[x,%{dx)]\x) ~dx{\ + -ix)\x), (8) 

where the error is in second order in dx. Then we get, 

= + 7x), (9) 

with the following uncertainty relation, 

n 



AxAp^ > -(l + 7(a;)), 



(10) 



which expresses the fact that the uncertainty in a mea- 
surement depends on 7 as well as on the average position 
of the particle. 

Next, it is straightforward to obtain an expression for 
the modified momentum operator in the cc-basis using the 
definition ([6|), 



(11) 



From (jlip the modified momentum operator can be writ- 
ten for short as p-y = —ifiD-y with 



Dy = (1+70;) 



dx 



(12) 



being a deformed derivative in space. In the 
^-representation, the corresponding time-dependent 
Schrodinger equation is 



d 

ih—{x\a,ts} = {x\H\a,ts), 
ot 



(13) 



and if we consider the Hamiltonian operator to be, H = 
pi^l2m+V{x), we arrive at the following differential equa- 
tion: 



^ 9 , , . 
tn-Q^V[x,t) = - 



2m 



D^^^P{x,t) + V{x)4>{x,t). (14) 



We now focus our attention on the case of a single 
spinless particle system. If the wave function 'ip(x,t) is 
normalized, it is possible to define a probability density 
p{x,t) = \4>{x,t)\'^. Using Ea. ((T4)) it is straightforward 
to derive a modified continuity equation 



where the probability flux is given by, 

Ti{l + 7x) f dip di(j 



(15) 



ax ax 



For standing waves in a null potential, the wave function 
(j){x) satisfying (fT4)) obeys 



2m 



DU{x) = Ecj>{x), 



(17) 



fi' d^ 



d ( I \ d(j){x) 



2me dx-^ 2 dx \2m£ 



dx 



+ E4){x) =0, (18) 



with me = m/(l -f 72^)^ being the particle's effec- 
tive mass, in perfect analogy with problems involving a 
position-dependent mass particle in semiconductor het- 
Here this particular expression for 



HI 



erostructures 

the effective mass arises naturally from the non-additive 
translation operator. Equation ()18p can then be rewrit- 
ten in the form of the Cauchy-Euler equation [l9| . 



2d'^4'{x) 



au{x)^^^+b<j>{x)^0. (19) 
du 



with u{x) ^ {1 + jx), a 1, and 5 2mE / {h^'f^). The 
general solution for Eq. p^ is 



(x) = exp 



±z— In (1 + 7a;) 
7 



(20) 



where fc is a continuous variable regarding the particle's 
wave vector. Although the wave function for a free par- 
ticle is not the usual plane wave, the particle's energy 
spectra is continuous, E = and independent of 7. 

For this free particle, the probability flux Jj is the same 
as an object moving at the classical velocity fik/m. 

If we now assume that the particle is confined to a 
one-dimensional infinite well of length L, the boundary 
conditions 0(0) — and (t>{L) = lead to the wave func- 
tion. 



An sin 
.0, 



■ In (1 -|- 72:) 



< a; < L 
otherwise. 



(21) 



where the wave vector is now quantized 



ln(H-7L) 



, with n= 1,2,3,4, ... (22) 



The energy for a confined particle can then be written 
as. 



En 



2mln^(H-7L)' 



(23) 



where ln(l+7L)/7 corresponds to an effective di- 
lated/contracted well that approaches L as 7 0. 

In Fig. [T] we show how the energy levels of a particle 
confined to an infinite well increase with n for different 
values of 7. Since the effective mass, me, in our de- 
scription, decreases with 7, and the lower the mass, the 
bigger is the kinetic energy of the particle, it means that 
the energy En increases with 7, as depicted in the in- 
set of Fig. [H The same reasoning can be applied to the 
size of the well, namely, the increase (decrease) of 7 above 
(below) zero, leads to a more pronounced contraction (di- 
lation) of its effective length. The asymmetry caused by 
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FIG. 1. The energy for a particle confined in an infinite well. 
The circles represent the energy for 7 = 0, the squares for 
7 = —0.5, and the triangles for 7 = 0.5. The energies are 
discrete, the solid lines are just guides to the eye. The inset 
shows the energy against 7 for the three lowest levels. 



the parameter 7 can be adequately quantified in terms of 
the average position of the particle in the box, calculated 
as (x) = (j)*x4>dx. From (PT|). we obtain the following 
expression: 



L 
2 



7' + 4fc2 



7 



(24) 



and the average of the modified momentum is {p^) = 0. 
Figure 2 shows the average position of the particle against 
7. As expected, when 7 = the average value is al- 
ways 0.5. The ground state is the most affected by the 
non-additivity of the space. As the quantum number 
increases, for example n = 20, the particle's average po- 
sition becomes independent of 7, (x) — ?> 0.5. 

The position operator in other space directions still 
commutes. Therefore, the theory developed here can be 
easily extended for two and three dimensions. For ex- 
ample, when considering a two-dimensional infinite well, 
the corresponding wave function can be expressed as the 
product y) = 4>i(x)4>2(y) , where 4)i(x) and 4'2(y) are 
the wave functions in the x and y directions, respectively. 
The contour plots for the probability density p{x, y) of a 
particle moving in a two-dimensional box are shown in 
Fig. Ela) rix ^ Hy ^ 1; in fig. ^h) for = 1, = 2 ; in 
Fig. ^c) for Tlx — Uy — 2; and in Fig. ^d) — Uy ~ 20, 
where we have used 7 = 1 in all panels. The ground state 
shows clearly that the particle spends more time out of 
the box center. These results also indicate that the cor- 
respondence principle remains valid, as for — riy — 20 
(large quantum numbers) the probability to find a parti- 
cle is practically the same everywhere in the square well 
[see Fig.Hd)]. 

From the above examples, for a null (free particle) or 



FIG. 2. The average position (x) of a particle confined in an 
infinite quantum well. The solid line gives the average value of 
X for n = 1, the dashed line is for n — 2, the case n = 3 is the 
dashed-dotted line, and the large quantum number n = 20 is 
represented by the dotted line. As n increases the asymmetry 
of the wave function is reduced, (a;) — 0.5. 



a constant potential (infinite ID and 2D well), the effect 
of 7 on the wave function of the particle is to stretch or 
contract the space variable. In this way, we expect the 
same effect for a particle subjected to a potential barrier 
with height Vb > 0, and located between x = and 
X = a. The wave function becomes a linear combination 
of Eq.®, 



+ re 



— ika 



for 2: < 



te 



ika' 



Be-^«°', for0<a;<a, (25) 
for X > a 



where k = ^^2mE/h, q = y2m{E - VQ)/h, and 
we used a' = In (1-1- 70)77 for short. The coefficients 
A, B,r,t are found taking the continuity of the wave func- 
tion and its spatial derivative in x = and x = a. The 
transmission and tunneling probability is modified by 7: 



1 



1 



sin^ qa' 
4E{E-Vo) ' 



4E{Vo-E) ' 



for E>Vo 



for E <Vn 



(26) 



Figure 4 shows the transmission probability against the 
energy ratio E/Vq for different values of 7. For E > 
Vb, we can see that the resonances (T = 1) depend on 
the value of 7. Increasing 7 is analogous to decrease 
the length of the barrier potential, and the transmission 
probability gets closer to 1. When E <Vo, the quantum 
tunneling probability increases with 7 (thinner barrier). 
The inset in Fig. 0] shows an oscillation with increase 
period for E > Vq. 



FIG. 3. The contour plot of the probabihty density for a 
particle in a two-dimensional box for 7 = 1, where the quan- 
tum numbers used are: a) Ux = riy — 1, b) = l,ny — 2, 
c)nx = Uy = 2 and d) nx = riy — 20. The probability in- 
creases from blue to red. 



FIG. 4. The transmission probability T for a rectangular 
barrier for three different values of 7. The solid line is for 
the normal space (7 = 0). The dotted curve is for 7 — 0.5, 
while 7 — —0.5 corresponds to the dashed curve. The inset 
shows the transmission probability for three different values 
of energy. Here we used ^/2mVo/h'^ ~ 6 



Let's turn our attention to the kinetic operator devel- 
oped through the non-additive approach introduced here. 
According to Eq.ljlip. we can write the modified momen- 
tum operator as pj — (1-f 7a;)p, so that the kinetic energy 
operator becomes. 



K = — 
2 . 



1 1 



(27) 



By comparison, this expression does not constitute a par- 
ticular case of the general kinetic energy operator pro- 
posed in Ref. 21| to describe a position-dependent mass 
in the effective mass theory of semiconductors, namely, 
I [m°'pm^pm^ + m^pm^pnf), with a + P + 6 — —1. To 
the best of our knowledge, and despite its claimed gen- 
erality in terms of the parameters a, /3, and 7, this last 
operator has not been deduced from first principle calcu- 
lations. 

In summary, we have introduced a non-additive trans- 
lation operator that can be identified as a g-exponential 
(1^ . By means of this operator, we have developed a 
modified momentum operator that naturally leads to a 
Schrodinger-like equation reminiscent of the wave equa- 
tion typically used to describe a particle with position- 
dependent mass. First, our results indicate that a free 
particle in this formalism has a continuum energy spec- 
trum with a wave function that is a modified plane wave. 
For a constant potential, like the problem of a particle 
confined to an infinite well, the energy now depends on 
the parameter 7, and for potential barrier, the peaks of 
maximum transmission probability are 7-dependent. In 
this context, we can argue that the substrate, as defined 
here, behaves like a graded crystal whose local properties 
determine the effective mass of the confined particle. Our 



approach can therefore be useful to describe the particle's 
behavior within interface regions of semiconductor het- 
erostructures. Another interesting consequence is that 
one can map this non-additive theory to an additive one 
using the appropriate potential. As future work, we in- 
tend to investigate the behavior of the non-additive par- 
ticle system when subjected to confining potentials that 
depend explicitly on position, for example, the harmonic 
oscillator or the central potential cases. 
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